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Abstract 

A cosmological model in which the primordial perturbations are 
provided by global monopoles and in which the dark matter is cold 
has several interesting features. The model is normalized by choosing 
its single parameter within the bounds obtained from gravitational 
wave constraints and by demanding coherent velocity flows of about 
600km/sec on scales of 50h~ 1 Mpc. Using this normalization, the model 
predicts the existence of dominant structures with mass 2 x 10 16 M Q on 
a scale 35h~ 1 Mpc i.e. larger than the horizon at t eq . The magnitude 
of the predicted mass function in the galactic mass range is in good 
agreement with the observed Schechter function. 



1 Introduction 

Recent observations have revealed the existence of nonlinear structures on 
very large scales (up to 50/i -1 Mpc). This observational fact has been a chal- 
lenge for theoretical models during the last few years. In fact the currently 
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popular models based on adiabatic primordial perturbations produced dur- 
ing inflation have difficulty to account for such large structures while main- 
taining successes on smaller scales. It is therefore important to investigate 
alternative theories. 

A new class of such theories based on primordial perturbations produced 
by seeds rather than a superposition of waves with random phases has at- 
tracted significant attention during the past few years. Physically motivated 
candidates for such seeds are topological defects jlj]. Topological defects are 
localized energy concentrations which are predicted by many particle physics 
models to form during phase transitions in the early universe. According 
to their geometry topological defects appear in the form of monopoles (sta- 
ble pointlike defects), cosmic strings (linear defects) domain walls (planar 
defects) and textures (unstable pointlike defects). The vast majority of lit- 
erature had until recently focused on gauged cosmic strings: linear defects 
interacting with short range forces. It was later realized[§] however that 
global defects interacting with long range forces, have two important advan- 
tages over gauged defects: 

1. They generically have a larger correlation length during their evolution 
and can therefore seed structures on larger scales than gauged defects. 

2. Their long range Goldstone field can provide efficient annihilation 
mechanisms that prevent them from dominating the energy density 
of the universe and lead to a scale invariant distribution: the scaling 
solution. 

Thus, for example, monopole-antimonopole annihilation makes global monopoles 
consistent with standard cosmology ||, in contrast to the case of gauged 
monopoles. These advantages have led to recent extensive study of the un- 
stable pointlike defect -global texture- with very encouraging results ||. 

In this letter we consider a model in which the large scale structure is 
seeded by stable pointlike defects (global monopoles) and in which the dark 
matter is cold. Global monopoles may form in particle physics models where 
the breaking of a global symmetry results in a vacuum manifold with topol- 
ogy S 2 . Such a breaking is realized in a model with the simple Langrangian 
density: 

C = ■ 8"$ - ~A(<& • $ - t} 2 ) 2 (1) 

where $ is a scalar field with 0(3) global symmetry and the symmetry 
breaking is 0(3) — * 0(2). Global monopoles are topologically nontrivial so- 
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lutions of the above model. They are described by the spherically symmetric 
ansatz: 



$ = /(r)- (2) 
r 

Even though a solution for f(r) is not known analytically, it is easy to 
show using the equations of motion that its asymptotic behaviour is f(r) — > 1 
as r — > oo and f(r) — > r for r — > 0. Thus the energy momentum tensor may 
be approximated at large distances by 

t/^t;^, If = 2#~0 (3) 

Solving the Einstein equations with this leads to the global monopole 
metric outside of the monopole core 0] 

ds 2 = dt 2 - cfr 2 - (1 - 8TTGrj 2 )r 2 {d9 2 + sin 2 0d0 2 ) (4) 

Thus there is a spherical deficit angle A = ^irGrj 2 . Equivalently any 
planar surface containing the monopole, has the geometry of a cone with 
deficit angle A. A more detailed analysis shows that in addition to the 
deficit angle the monopole induces an effective small negative mass M in its 
spacetime [11]. The magnitude of this mass however is proportional to the 



scale of symmetry breaking rj and is too small to be relevant for structure 
formation. On the other hand, the effects of the deficit angle are independent 
from the distance to the monopole core and may therefore be important even 
for very small A. By obtaining the geodesic equations in the spacetime (4) 
it may be shown that a monopole moving with velocity v m induces a velocity 
perturbation Av towards the line of its trajectory, to its surrounding matter. 
Ignoring terms of 0( Vmatter ) and 0( " matter ), the magnitude of the dominant 

^monopole C 

component of Av is Av = <iTTGr] 2 v m ^ m and its direction is perpendicular 
to the monopole trajectory . These primordial velocity perturbations may 
then grow to form galaxies, clusters and large scale structure. In what 
follows we obtain some of the predictions of the above described model. We 
make the assumption that the monopoles move with relativistic velocities 
on straight line trajectories. This assumption is justified due to the long 
range attractive monopole- antimonopole forces ||],@ which are expected to 
induce relativistic velocities. 

The structure of this paper is as follows: In the next section we use the 
Zeldovich approximation to obtain the growth of perturbations induced by 
global monopoles. Then we obtain the peculiar velocities predicted by the 
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model and we use this result to normalize the single free parameter rj of 
the model. In section 3 we use the results of the Zeldovich approximation 
to obtain the spectrum of nonlinear masses predicted by the model. We also 
obtain the predicted mass function and we compare the galactic mass range 
with the observed Schechter function. Throughout the paper h is the Hubble 
constant in units of WOkm/ (sec- Mpc) and we consider a spatially flat = 1 
universe. 



2 The Zeldovich Approximation 

We will use the Zeldovich approximation to calculate the growth of the 
velocity perturbations induced by relativistic global monopoles. The cal- 
culation is similar to the one in the case of cosmic string wakes [12|, ]|] 
with the important difference that in this case the geometry of each wake 
is cylindrical rather than planar. Consider a particle with unperturbed co- 
moving distance q from the trajectory of a global monopole (assumed to 
be a straight line). Due to the initial velocity perturbation induced by the 
monopole a comoving position displacement ^(q,t) will develop. Thus the 
physical distance r±_ from the monopole trajectory may be written 

r ± (q,t) =a(t)(q-*(q,t)) (5) 

where a(t) is the scale factor of the universe. Also, in the Newtonian ap- 
proximation we have 

r± = - Vr ± $ (6) 
V^<*> = 4:irGp (7) 

Using (5), (6), (7) and linearizing in ^ leads to 

^ + 2-^ + 3-^ = (8) 

a a 

The velocity perturbation induced by a monopole at time ij is Av = 
4irGv m j m . Therefore, the initial conditions to be used for the solution of 
(8) may be written in the comoving frame as: 

= 0, *(* i ) = Au(^)f J U>t eq (9) 

= 0, *&) = AiA§(^)3, t t <t eq (10) 
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The evolved comoving displcement ^f(t,ti) with t > t eq is now easily 
obtained from (8), (9) and (10) and keeping only the growing mode it may 
be written as 

3 ^ ^ 2 

¥>(t,ti) = -Au^r'ts, (11) 
o 

*<(t,ti) = ^A^4tft^e7(l + ln(^)), t 4 <t eg (12) 

This result may be used to calculate several interesting quantities: 

1. The large scale peculiar velocities predicted by the model. 

2. The nonlinear mass of objects seeded at an initial time ij. 

3. The mass function (for this we will also need the scaling solution for 
monopoles). 

Here we find the predicted peculiar velocities and use this result to normalize 
the single free parameter of the model r\. 

It is easy to show using (11) that the predicted peculiar velocity caused 
by a single monopole, coherent on a scale q is 

v>(t ) = tf>(t ) = r(Au)(^)S = — (Gr] 2 ) 6 v m j m q^km/sec (13) 

« H 1/3 

where q$o is the comoving scale in units of 50/i _1 Mpc, u > lj is the number 
of monopoles per horizon volume in the scaling solution (tf om = vzq) and 
(Gr? 2 ) 6 is Gr? 2 in units of 10~ 6 . It may be shown using arguments similar to 
those in Ref.|| that the effect of all later monopoles is to produce an rms 
velocity larger by a factor of 1.64 than the result (13). Typically v3 ~ l||. 
Demanding v > (to)\ q50 ~ 600km/ sec we find (Grj 2 )QV m ^ rn ~ 10 implying, for 
Vmjm — la symmetry breaking scale rj for the monopole forming phase tran- 
sition rj ~ 3 x 10 1G GeV. This value of 77 is consistent with the gravitational 
wave background constraints obtained in Ref.[|ljJ rj < 5 x 10 ie GeV. 

3 Nonlinear Structures 

The solution (11) , (12) describes how the expansion of r±(t) is slowed down 
gravitationally due to the monopole induced velocity perturbation. The 
Zeldovich approximation used in deriving it can be used until the shell de- 
scribed by r±(t) turns around i.e. r±(q,t) = The condition r±(q n i,t) = 
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determines the thickness q n i of the cylindrical overdensity formed by the 
initial velocity perturbation. Using (5) it is easy to show that q n i(t,U) = 
2ty(t,ti). The nonlinear scale q n \ however can not grow beyond the comov- 
ing scale of the initial velocity perturbations tf m . Therefore, the correct 
expression for the turnaround scale at the present time to is Qnl{to,U) = 
min(^ om ,2^(t ,ti)) or, defining t c ° m = 2^(t ,Q we have: 

qnl(t ,U) = C m , h<t* (14) 
q n l(to,U) = 2^(t ,ti), U>U (15) 

The form of q n i(to,ti) is shown in Figure 1 for > t eq . Notice that for 
ti < t* it is t? om , not 2^(to, U), that determines the scale q n \. It is straightfor- 
ward to calculate t* and show that = t eq f(X) where A = 0A(G / q 2 )6h 2 v r n'jm 
and 



/(A) = 6^, A<1 
/(A) = A§, A>1 

In what follows we will use the normalization of the previous section and 
assume A = 4/i 2 > 1 

The nonlinear mass corresponding to structure originating at time ti may 
be written as 

M(t, U) = 7rq 2 nl (t, ^> m C m p(*o) (16) 
For A > 1 we have t* > t eq and using (14), (15) and (16) we obtain 

M(t ,ti) = U<t eq (17) 

M(t ,U) = ^U, t eq <t i <U (18) 
M(t ,U) = |^^(|)^o, U>U (19) 

For ti < t*, q n i(to,ti) ~ i? om and therefore all matter goes nonlinear on 
structures seeded before i*. However, M(io, U) is an increasing function of ti 
and therefore structures seeded at ti < i* can form by accreting structures 
seeded before ti. On the other hand, for ti > t*, M(to,ti) decreases with 
ti and since all matter has gone nonlinear in structures seeded before i*, 
monopoles after t* can not produce nonlinear structures by today but can 
only induce large scale velocity flows. 
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The important point to notice in Figure 1 is that the distinguished scale 
on which the largest and most prominent structures form is t* i.e. a scale 
which is larger than ig° m by a factor A5. Thus the characteristic scale of 
dominant structures in the model is 



tf m = \hl° q m = Zhh- x Mpc (20) 

while the corresponding mass is 

M(to,U) ~ 2 x 10 16 /i _1 M Q (21) 

where we have used the normalization obtained in the previous section A = 
Ah 2 . 

It is straightforward to now obtain the mass function predicted in the 
model using the additional input of the monopole scaling solution. Since 
the number of monopoles per Hubble volume per conformal horizon n is 

= ^ with v ~ 0(1), it is easy to see that the number density of nonlinear 

structures with mass larger than M is given by 

n>M = / dM^ff- (22) 
Jm t?{M) 

where M(io>ij) was obtained in (17)-(19). It is straightforward to compute 
the integral (22) to find n>M- The result is shown in Figure 2 and may be 
written as 



We also show the Schechter mass function (obtained from the Shechter lu- 
minocity function|i~0|] with a mass to light ratio of 100 ; ^^) in Figure 2 for 
comparison with the monopole result. Clearly the model predicts a cutoff 
at the mass scale M* . This cutoff however is much larger than the galactic 
mass and should be distinguished from the observed cutoff in the Schechter 
function. The Schechter function cuttof is expected to naturally emerge 
when effects of hydrodynamics are taken into account as was the case in 



Ref. HJ. In particular [14], cooling effects may provide the cutoff since the 
gas on larger scales can not cool in one Hubble time. As shown in Figure 2 
the magnitude of the predicted mass function in the galactic mass range is in 
reasonably good agreement with the Schechter function given the assump- 
tions involved in the calculation. The slope however is —1 for the monopole 
model but only 0.3 ±0.1 for the observed Schechter function. This potential 
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problem of the model could be resolved by considering hot dark matter or 
by introducing bias in the model. 

We have therefore presented a study of a model for large scale structure 
formation based on primordial perturbations created by global monopoles 
and nonrelativistic dark matter (cold dark matter). The model has several 
encouraging features and clearly deserves further study. The model's single 
free parameter was fixed by demanding agreement with the observed large 
scale velocity flows. The obtained value was consistent with gravitational 
wave background bounds. The characteristic scale of dominant structures 
in the model is larger than the comoving horizon at t eq implying that the 
model has significant power on large scales. The main distinct feature of 
the monopole wakes considered here when compared to the cosmic string 
wakes is their generic cylindrical geometry to be compared with the planar 
geometry of string wakes. In addition, the global nature of the field in 
the monopole model implies the existence of larger correlation length in the 
scaling solution compared to the case of cosmic strings. Thus the parameters 
of the scaling solution in the case of monopoles favor the formation larger 
structures than cosmic strings. 

One of the potential problems of the model is its prediction that all 
mass has gone nonlinear at the present time to- This problem, which also 
appears in the case of cosmic strings with cold dark matter could be resolved 
by considering relativistic particles as the dark matter (hot dark matter). 
Work in this direction is in progress. The consideration of hot dark matter 
could also soften the slope of the mass function which is predicted to be 
steeper than the Schechter function for galaxies. 

The global monopole model for structure formation predicts a distinct 
signature on the microwave backgroud. In particular, a global monopole 
moving perpendicular to the line of sight of the observer will produce a tem- 
perature dipole in the microwave sky i.e. a hot-cold spot pair. The detailed 
form, magnitude and distribution of this signature will be presented in a 
separate publication 
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Figure Captions 

Figure 1: The nonlinear mass q n \ vs time of initial perturbation tj. 
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Figure 2: The predicted mass function superposed with the observed Schechter 
function (obtained by using a fixed mass to light ratio). 
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